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RESUMEN 

La Termodinámica Química es la teoría unificada de las transformaciones físicas y químicas. Este 
trabajo presenta tres teoremas sobre los cuales se puede fundar toda la disciplina. El primero de 
ellos revela la conexión existente entre el método de los multiplicadores de Lagrange de la 
formulación no-estequiométrica y la condición de homogeneidad implícita en la mayoría de las 
ecuaciones termodinámicas. Sin embargo, esta última es menos restrictiva que la condición de 
homogeneidad del sistema, es decir, de todas sus propiedades; se refiere exclusivamente a la 
homogeneidad del potencial de Gibbs respecto del vector composición. Por ello, las ecuaciones 
termodinámicas también pueden ser aplicadas a los sistemas multifásicos en equilibrio. El 
segundo indaga en la definición de elemento empleada en la aproximación física a la cuestión del 
equilibrio químico: los potenciales de los elementos no dependen del estado de las partículas 
consideradas porque el potencial de Gibbs es una función homogénea del vector de las 
composiciones de los elementos. El tercero muestra la relación existente entre la definición del 
potencial de los elementos y la dinámica del equilibrio a nivel microscópico: esta definición 
corresponde al cambio en el número de partículas en cualquier estado — los potenciales de los 
elementos son magnitudes homogéneas de grado cero — . Los tres teoremas son maneras 
diferentes de presentar la misma idea: el concepto de homogeneidad y sus consecuencias. 

Palabras clave: potencial de Gibbs; Termodinámica Química; multiplicadores de Lagrange; 
transformaciones de Legendre; equilibrio dinámico 
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1. INTRODUCCIÓN 

El éxito de Gibbs al resolver el problema de la energética de los cambios de fase 
llevó a transferir sus conclusiones a las transformaciones químicas, lo cual dio 
origen a la Termodinámica Química. 

En las transformaciones físicas, las dos formas (estados) de una sustancia están 
separadas por una superficie. En las transformaciones químicas, ambas formas 
están mezcladas. La única diferencia entre el cambio de fase de un cuerpo puro y 
una isomerización en un medio homogéneo está en la energía libre de mezcla. En 
el caso de una transformación líquido-vapor, el equilibrio es posible porque la 
energía de las partículas — que es mayormente cinética en el estado de vapor — es 
mayormente potencial en el estado líquido. En el caso de la isomerización, la 
coexistencia de ambas formas se debe a la estabilización que produce la mezcla, a 
la cual siempre está asociada una energía libre negativa. 

Estas diferencias no impiden el desarrollo de una formalización común para todas 
las transformaciones, lo cual es posible al trabajar con propiedades molares 
parciales: el efecto de la mezcla correspondiente a un cambio diferencial de la 
composición es nulo. Por lo tanto, la condición de igualdad del cambio en el 
potencial de Gibbs asociado a un cambio de la masa del elemento (¡j) es la 
condición de equilibrio para cualquier transformación, física o química. 

En este trabajo se indaga en las razones profundas para la unificación del estudio 
de todas las transformaciones. Con este propósito, se demuestra que la condición 
de homogeneidad del potencial de Gibbs respecto del vector composición es la 
responsable de la existencia de los potenciales de los elementos (multiplicadores) 
(Primer Teorema). La existencia de dichos potenciales hace caer la barrera que 
separa a las transformaciones físicas de las químicas, es decir el cambio de 
nombres que acompaña a las últimas. Sobre la base de esta conclusión, es posible 
asignar un cambio de nombres a las transformaciones líquido-vapor o no hacerlo 
en el caso de las isomerizaciones (Segundo Teorema). La igualdad de las 
propiedades molares parciales ¡u lleva implícita la idea de dinámica a nivel 
microscópico, lo cual se demuestra por medio de experimentos (Tercer Teorema). 

El Primer Teorema revela la conexión existente entre la tradicional formulación 
no-estequiométrica de los problemas de equilibrio químico (Smith y Missen 1991) 
y una variante del método introducido por Alberty (1997) para obtener la 
expresión de los potenciales de los elementos. El Segundo Teorema muestra que 
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los potenciales son independientes del estado, y es una generalización del 
enunciado de Beattie y Oppenheim (1979). El tercer teorema resalta el aspecto 
dinámico de la definición de potencial de un elemento (Gerasimov 1971). 

Los tres enunciados están íntimamente relacionados — por lo cual se los presenta 
juntos y en ese orden — y pueden servir de base para fundar toda la 
Termodinámica Química (Luetich 2001). 

Este artículo es el tercero de un grupo de cuatro trabajos relacionados escritos 
para mostrar el sentido físico del concepto de reversibilidad microscópica. Aquí, 
el énfasis se pone en los conceptos de homogeneidad y dinámica a nivel 
microscópico. 

2. ENUNCIADOS Y DISCUSIONES 

PRIMER TEOREMA 

MINIMIZACIÓN LINEAL Y LA CONDICIÓN DE HOMOGENEIDAD 

Cuando la función objetivo del método de los multiplicadores de Lagrange es una 
función homogénea de las variables sujetas a restricciones, la función auxiliar 
(lagrangiano) es la transformada de Legendre completa de la función objetivo 
respecto de las mismas variables. La función auxiliar y la función objetivo 
contienen la misma información. 

HIPÓTESIS 

a. Sea la función 

G = G(m, n) (2.1) 

homogénea respecto de n, es decir, 

G(m, kn) = kG(m, n) : keR. (2.2) 

Además, 

c 

G = Yj n ^i' < 2 - 3 ) 

¿=i 

siendo n¿ las componentes de n, y C la dimensión de n. 

b. Existe una matriz C = CexC , tal que 
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n* = YjCijm, j = 1,...,E<C, (2.4) 

i=\ 

donde c¿,- es la componente de C correspondiente a la fila i y a la columna j, y E es 
la dimensión de n* . 

DEMOSTRACIÓN 

De acuerdo a las hipótesis a y b, G podría ser expresado en una forma más 

compacta: 

E 

G = J]n*{ j , (2.5) 



donde los £,- son combinaciones lineales de los ju¿ . 
Definiendo la función transformada 

E 

(2.6) 



G* ± G-^nHj, 




se tiene: 




C E 

Z ^ ~ Z n *j & 
í=i j=i 


= G* 



(2.7) 



Entonces, sustituyendo las variables n*, se obtiene: 

C E C 

Z^ ¿ -ZZ c ^ ¿ £ = G *- (2 - 8) 

¿=1 j=\ 1=1 



De una forma más práctica, 

rit = G*, (2.9) 



¿=1 ^ ^ ' 



y, dando un nombre al primer factor de la sumatoria, 

c 



Z^ ¿ = G*. (2.10) 

¿=i 

Para que esta igualdad sea válida cuando G* — 0, 
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ti = o, 

y, finalmente, se deduce que 

E 

fr-J]cij£j= O, i = 1, ..., C, 

expresiones de las cuales se pueden obtener los £¡ . 
Paralelamente, de acuerdo a la hipótesis a, 

— ñC — 

dG = ^j -z dm k + ¿] # drii, 



k=i 



dm k 



(2. 11) 



(2. 12) 



(2. 13) 



i=i 



donde S es la dimensión de m, y 

c 
(dG) m = ^jUidrii. 
¿=i 

De la ecuación (2. 6), 
(dG*) m = 



E C 

dG - Yj Yj & d (c y Ui) + Ci J n¿ d ^ 
j=i ¿=i 



dG* \ 

drii ) 



= fii 



E C 



mn[¿] 



j=i ¿=i 



6 



dtt,; 



+ C í/ -Al ¡ — ^ 



mn[¿] 



i = 1, ..., c. 



(2. 14) 



(2. 15) 



(2. 16) 



esta ecuación es equivalente a la ecuación (2. 12), dado que el primer miembro es 
cero y el segundo término de la suma es nulo. 



Cuando G* = 0, 

(dG*) m = 0. 
Por lo tanto, 



(2. 17) 



E 



(dG*) m = -J]n*d^ 

7=1 



(2. 18) 



De este modo, por un lado, de acuerdo a las ecuaciones (2. 5) y (2. 18) y a la 
definición (2. 6), G* es la transformada de Legendre completa de G con respecto a 
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las componentes de n*; y, por otro lado, de acuerdo a la misma definición, G* 
puede ser vista como la función auxiliar de Lagrange (lagrangiano) para la 
optimización de G sujeta a las condiciones (2. 4): £,- son los "multiplicadores". 

COMENTARIOS 

En el método de los multiplicadores de Lagrange, las variaciones de G son 
combinaciones lineales de las variaciones de n. Ésta es la misma razón por la 
cual decimos que G es una función homogénea de primer orden de n (figura 1): la 
linealidad y la homogeneidad son propiedades equivalentes. 
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Figura 1 

Ilustración de la homogeneidad de G respecto de n en un sistema binario líquido-vapor. 

Para la región encerrada: (i) G(ni) y (ii) G(Ani) = AG(ni). (Aquí, A = 2.) 

(G no es homogénea respecto de V) 

Por lo tanto, cuando se realiza la transformación de Legendre completa de una 
función a la cual se aplican las conclusiones del teorema de Euler (función 
homogénea), los resultados coinciden con los que se obtienen usando los 
multiplicadores de Lagrange. 
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La función auxiliar de Lagrange puede ser identificada con la transformada de 
Legendre completa de la función G con respecto a las componentes de n* 
(Luetich 2002b). Nótese que el lagrangiano 



$ = fc(n,¿), 



(2. 19) 



y la transformada de Legendre completa de G(m, n) con respecto a n es: 

G* = G* (m, n), (2. 20) 

pues m permanece constante y la transformación de Legendre se hace con 
respecto a n* , 



G* = G*(m,¿) = G*(0 = G*(n*,Q. 



(2.21) 



Los pares de funciones {G(m, n), G*(m, ¡x)} y {G*(m, 0, G(m, n*)} contienen la 
misma información porque cada una de ellas es la transformada de Legendre 
completa de la otra respecto de la segunda variable. Para relacionar estos dos 
conjuntos y mostrar que G(m, n) contiene la misma información que G(m, n*), 
este teorema muestra que G* (m, ¿) contiene la misma información que G(m, n). 



transformación de Legendre 
completa respecto de n 




transformación de Legendre 
completa respecto de |j. 




n* = A n 



G-(n*) T ^ 



n* = An 



1 
1 


transformación de Legendre^\^ 
completa respecto de £ 


i 
i 


G(m, n*) 


G*(m, Q 




transformación de Legendre 
completa respecto de n* 







Figura 2 

Demostración diagonal de la equivalencia de G(m, n) y G(m, n* ) 

Esta prueba "diagonal" (figura 2) muestra la equivalencia de la formulación 
no-estequiométrica usando multiplicadores de Lagrange y la formalización de la 
Termodinámica Química recurriendo a transformaciones de Legendre. 
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SEGUNDO TEOREMA 

LAS VARIABLES DE COMPOSICIÓN Y LOS POTENCIALES DE LOS ELEMENTOS 



En los sistemas en equilibrio, el potencial calculado para una partícula es 
independiente de su entorno. 

HIPÓTESIS 

Sea un sistema en equilibrio térmico, mecánico y químico; y sea a cierto estado 

de agregación de las partículas del elemento j-ésimo. 

DEMOSTRACIÓN 

Volviendo a la ecuación (2. 5) y tomando en cuenta que las variables 

transformadas incluyen a las partículas en sus varios estados, se obtiene 



n j = n 



7 + E*?' 



(2. 22) 



separando la parte correspondiente al estado a de las demás (<p) (figura 3). 




Zcp / 




\ y 



Ch 



$ 












Figura 3 

a y los otros estados (Z<p) del elemento A(#) en el sistema homogéneo de dos elementos 

A-B-AB-A 2 B 

El potencial de Gibbs de cada especie es la suma de los potenciales de Gibbs de la 
especie en sus distintos estados. 



G¡ = G" + Y G? 



"j ~ ^j 






(2. 23) 



El potencial de Gibbs de todo el sistema es la suma de los potenciales de Gibbs de 
los distintos estados. 
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E 

G = J^Gj = G a + Y ¡ G ' f ( 2 - 24 > 

7=1 f 

La variación del potencial de Gibbs con el número de partículas del elemento j en 
el estado de agregación a es: 



8G\ (8G a 

+ 



dn 1 'tph^ ] ydn] ) TPl ¿m 



( d^G*^ 

dn°¡ 



9G ° ■ ,2.25) 



TPn 



e\j\ 



dn< j 'TPn^ 



'<■[./! 



donde n^ representa al conjunto de las variables de composición, excepto la 
correspondiente al estado a del elemento j. 

El supraíndice a puede representar cualquier estado de agregación, es decir, 

fMm jm -if.) = „. 

V arij J TPll ¡m V orij J TPn m \ arij J T pn[j] 
Esta propiedad conduce a la definición del potencial del elemento j-ésimo: 

donde n*[j] representa a todas las componentes de n* excepto laj-ésima. 

Este potencial es el multiplicador de Lagrange del elemento de acuerdo a la 
transformación de las variables de composición (2.22), es decir la variable 
conjugada de la variable de composición resultante. 



El potencial por partícula es 
Cj _ ( dG 



J- {dN* j/Tmij] 



(2. 28) 



donde d- es el número de Avogadro. Su valor es el mismo para todas las 
partículas del elemento j. 

COMENTARIOS 

a y, en general, (p son los estados (de agregación) de las partículas del elemento 
j-ésimo. Los distintos estados pueden ser de agregación física o química. Por lo 
tanto, las transformaciones se representan usando los mismos símbolos en ambos 
miembros. Las especies que no cambian durante las transformaciones se llaman 
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elementos. Los ejemplos presentados en la tabla indican qué es lo que se debe 
entender por elemento y estado (Luetich 2002a). 

Table 

Elementos y estados en cinco sistemas: (i) heterogéneos, (ii)-(iv) homogéneos 





sistema 


estados de A 


estados de B 


estados de C 


(i) 


A(h'q) ^ A(vap) 


líquido | vapor 








(ü) 


2A^A 2 


libre | A 2 








(iü) 


A + B ^ AB 


libre | AB 


libre | AB 





(iv) 


r A + B ^ AB 
l2AB^B + A 2 B 


libre | AB | A 2 B 


libre | AB | A 2 B 





(v) 


AC + B ^ C + AB 


AB|AC 


libre | AB 


libre | AC 



TERCER TEOREMA 

LA DINÁMICA DEL EQUILIBRIO Y LA DEFINICIÓN DE LOS POTENCIALES DE 

LOS ELEMENTOS 

En los sistemas que no muestran cambios a nivel macroscópico, los cambios 
posibles a nivel microscópico son aquellos que producen variaciones iguales y 
opuestos en los potenciales de Gibbs correspondientes. 

HIPÓTESIS 

Sea un sistema en equilibrio mecánico, térmico y químico. 

DEMOSTRACIÓN 

Supóngase que una cantidad no-observable de átomos del elemento j-ésimo 
(partícula) pasa del estado de agregación a al estado de agregación /? (figura 4). A 
esta transformación estaría asociado el siguiente cambio en el potencial de Gibbs: 

(dG^ 



(dG) 



TPn' 



í 8G a 

V dn j 'TPn 



dn a + 



\ dn) , 



dn p . 



(2. 29) 



TPn' 



Debido a que el potencial de Gibbs del sistema puede ser expresado como sigue: 

G = G a + G p + Yj GiP > ( 2 - 3 °) 



y que los estados <p no están involucrados, 
dG = dG a + dG p . 



(2.31) 
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Figura 4 

Dinámica a nivel microscópico. Variaciones opuestas de la masa del elemento j-ésimo 

en dos estados (a y J3) en equilibrio 



Recordando el resultado del Segundo Teorema y la propiedad de las derivadas 

dG\ í -dG\ ( dG 



dn 
se obtiene 

dG a 



- dn 



-dn 



(2. 32) 



d n J J TPn* 



dCr B \ 

\ dn] , 



(2. 33) 



TPn* 



Dada la relación existente entre los conceptos de potencial y fugacidad, esta 
ecuación es la expresión matemática de la definición de fugacidad como la 
capacidad de los "átomos" de escapar de un cierto estado. 



Por lo tanto, si 



dn] = - dn a j, 



se puede obtener la relación 

8G a = - dG 13 



(2. 34) 



(2. 35) 



que constituye la formulación termodinámica de la idea del equilibrio dinámico. 

Esta propiedad es una consecuencia de la igualdad de los multiplicadores de 
Lagrange de las partículas de un mismo elemento en dos estados distintos. 

( dG 



^ ] (dN° 



d = ^ 



TPN* 



K dN¡, 



(2. 36) 



TPN* 
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El efecto del movimiento de una partícula no observable sobre el potencial de 
Gibbs de todo el sistema es nulo. 

(dG) TPn . = (2.37) 

COMENTARIOS 

En el Primer Teorema se demostró que la equivalencia entre la linealidad y la 
homogeneidad da origen a la igualdad de los potenciales (propiedades 
homogéneas de orden cero). 

El orden cero de la homogeneidad de los potenciales de los elementos ¡j.: 

fi (x, kn) = ¡j. (x, n), (2. 38) 

donde x es el vector de las fracciones molares n/Z¿ n¿ , se debe al hecho de que el 
potencial del elemento es la derivada con repecto a una variable de una función 
homogénea de primer orden con respecto a esa misma variable. 

G(x,kn) = kG(-x,n) (2.39) 

dG(x,kn) d[kG(-x,n)] kdG(x,n) <9G(x, n) 



dkn kdn kdn dn 



(2. 40) 



En un sistema homogéneo, las dos variables son homogéneas, pero a la palabra 
homogéneo se le debe asignar un sentido distinto: G no cambia con el número de 
partículas que se consideren (o sea, no cambia con las unidades: a un mol 
corresponden tantas calorías como kilocalorías a un kilomol); m, en cambio, tiene 
el mismo valor en todos los puntos. 

A esta altura, es importante destacar que las propiedades molares parciales 
tienen el mismo valor que las propiedades molares no sólo en el caso de los 
cuerpos puros sino también cuando el tamaño del sistema es tan grande que el 
agregado o la quita de una pequeña cantidad de un componente no altera la 
composición. 

En el nivel atómico-molecular, la notación de Leibnitz (y todo el cálculo 
infinitesimal) resultan especialmente reveladores. 

En resumen, en un sistema en equilibrio, el movimiento de partículas de un 
estado de agregación a otro no está en desacuerdo con las leyes termodinámicas 
de los sistemas macroscópicos. 
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La idea de la dinámica del equilibrio está estrechamente relacionada con el 
concepto de reversibilidad microscópica (Luetich 2002c). 

3. CONCLUSIONES 

La condición de homogeneidad de G con respecto a n es el punto de partida del 
razonamiento que sucesivamente conduce a: (1) la definición de potencial de un 
elemento; (2) la aproximación física al problema del equilibrio; y (3) la 
formalización de la dinámica del equilibrio a nivel microscópico. Esto es así 
porque: (1) la homogeneidad está relacionada con el carácter lineal del método de 
los multiplicadores de Lagrange; (2) los multiplicadores de Lagrange son las 
variables conjugadas de las variables de composición de los elementos; y (3) los 
cambios en el estado de agregación a nivel microscópico no modifican el potencial 
de Gibbs del sistema porque, en el equilibrio, los potenciales de los elementos son 
propiedades homogéneas de grado cero. 

La condición de homogeneidad no es tan restrictiva como se podría pensar 
considerando el hecho de que llamamos homogéneo a un sistema cuando sus 
propiedades son todas homogéneas (fase). Si tomamos cualquier sistema simple 
en equilibrio y consideramos un número de partículas rc¿ de cada clase, el 
potencial de Gibbs de ese subsistema es G; si consideramos una cantidad Arc¿, el 
correspondiente potencial de Gibbs es XG, independientemente del estado de 
agregación de las partículas (figura 1). Por lo tanto, la misma razón que lleva a la 
definición de los potenciales de los elementos, es decir, la homogeneidad de G con 
respecto a n, también da sustento a la aproximación física a los problemas de 
equilibrio. 

En los sistemas con muchos elementos, el potencial de cada uno de ellos es la 
derivada del potencial de Gibbs con respecto a la variable de composición 
correspondiente, manteniendo las restantes variables constantes. En términos de 
variaciones (notación de Leibnitz), el potencial de un elemento es el cociente que 
resulta de dividir el cambio del potencial de Gibbs por el cambio de la variable de 
composición: la idea de variación, es decir de dinámica — aun en sistemas 
cerrados — está implícita en la definición de potencial. Por ello, las variables 
restantes deben permanecer constantes y los efectos del mezclado impedidos. 
Pero precisamente las variaciones diferenciales son las que no producen estos 
efectos, debido a que no modifican la composición global. No sólo la dinámica a 
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nivel microscópico no está prohibida por la formaliz ación de la Termodinámica 
Química: está implícita en ella. 
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